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We theoretically study dynamic scalar polarizabilities of the ground and select long-lived excited 
states of dysprosium, a highly magnetic atom recently laser cooled and trapped. We demonstrate 
that there are a set of magic wavelengths of the unpolarized lattice laser field for each pair of states 
which includes the ground state and one of these excited states. At these wavelengths, the energy 
shift due to laser field is the same for both states, which can be useful for resolved sideband cooling 
on narrow transitions and precision spectroscopy. We present an analytical formula which, near 
resonances, allows for the determination of approximate values of the magic wavelengths without 
calculating the dynamic polarizabilities of the excited states. 

PACS numbers: 31.15.am, 32.70.Cs, 31.30.jg, 37.10.Dc 
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I. INTRODUCTION 

The dysprosium atom has many unique features which 
makes it useful for studying fundamental problems of 
modern physics. This is a heavy atom which has many 
stable Bose and Fermi isotopes (from A = 156 to A = 
164) and a pair of almost degenerate states of opposite 
parity at E=19798 cm -1 . These features were used in 
study of the parity non-conservation (PNC) [lJ-IH and 
possible time-variation of the fine structure constant @- 

Fermionic Dy has the largest magnetic moment among 
all atoms, and only Tb is as magnetic as bosonic Dy. This 
opens important opportunities in studying strongly cor- 
related matter when gases of Dy atoms is cooled to ultra- 
cold temperatures [13[ . Recent progress in Doppler and 
sub-Doppler cooling is an important step in this direc- 
tion [L3l - [l7| . In addition to narrow- line magneto-optical 
trapping (MOT) (l8| . further cooling on narrow opti- 
cal transitions might be possible using resolved-sideband 
cooling [H [20|. 

In this method, vibrational states of the atom may cou- 
pled such that successive photon absorption and sponta- 
neous emission cycles reduce the vibrational quanta by 
one, until the atoms are in the motional ground state 
of their optical potential [l9| . It is important that this 
resolved-sideband cooling is performed at magic wave- 
length of the laser lattice field [2l|, [22j. At this wave- 
length, the energy (AC Stark) shift due to laser field 
is the same for both states used in the cooling. This 
results in a trap potential that is the same for both 
states, and optical transitions between vibrational states 
can be well resolved. This allows spectral selection of 
cooling transitions — those which remove one vibrational 
quanta — without contamination by heating transitions 
which add vibrational quanta. Other benefits to optical 
trapping at magic wavelengths include enhanced preci- 
sion spectroscopy and longer-lived quantum memory for 
quantum information processing (QIP) [2l| . 



In this paper we calculate dynamic polarizabilities of 
the ground and three long-lived excited states of Dy and 
present a number of magic wavelengths for the transitions 
between them. We also present an analytical formula 
which allows the determination of approximate values of 
the magic wavelengths near resonances without calculat- 
ing the dynamic polarizabilities of excited states. The 
optical field is assumed to be unpolarized, though we es- 
timate that polarization would induce only small shifts 
in the magic wavelengths. 



II. CALCULATIONS 
A. Ab initio calculations 

The dynamic scalar polarizability a a of atomic state a 
is given by (we use atomic units: % = 1, m e = 1, |e| = 1) 



3(2J a + 1) 
1 

E a — E n — uj 



1 



E a - E n + UJ 

a\\B\\n) 2 , (1) 



where J a is total momentum of state a, E a is its energy, 
D = — J2i r i is the electric dipole operator. Summation 
goes over complete set of excited states n. 

We use the relativistic configuration interaction (CI) 
technique described in our previous papers 0, to 
perform the calculations. The single-electron and many- 
electron basis sets, the fitting parameters and other de- 
tails of present calculations are exactly the same as in 
Rcf. 5]. This simple method provides a good accuracy 
for low lying states of a many-electron atom. However, 
it does not allow for the saturation of the summation 
in Eq. ((T|) over a complete set of many-electron states. 
On the other hand, the contribution of the higher-lying 
states in the dynamic polarizability does not depend on 
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frequency at small frequencies. Therefore, for small fre- 
quencies we can rewrite Eq. (JT|) as 



3(2J a 
1 



— Y 



i 



E„ — En 



E n , — E„l 



(a||D||n'> 2 



(2) 



where the summation is over a limited number of low- 
lying near-resonant states and a constant a a is chosen in 
such a way that Eq. @ at uj = provides the correct 
value of the polarizability. 

Dysprosium ground state static polarizability is known 
to be 166 a? B [2J|. Static polarizabilities of excited states 
are not known and need to be calculated. We use an 
approximate approach in which the dysprosium atom is 
treated as a closed-shell system and the effect of electron 
vacancies in the open shells is taken into account via 
fractional occupation numbers. The static polarizability 
of a closed-shell system is given by 



2 ^ (c||D||n) 2 

O t c fc n 

cn 



(3) 



where the summation is over a complete set of single- 
electron states including states in the core (c) and states 
above the core (n). Electric dipole matrix elements are 
calculated using relativistic Hartrec-Fock and Hartree- 
Fock in external field approximations (25|. Note that 
core polarization needs to be included only in one of two 
electric dipole matrix elements in (|3]) (see, e.g., Ref. [26| 
for details). 

We use the standard B-spline technique [27j to gener- 
ate a complete set of single-electron states. An additional 
term is included into the Hartree-Fock Hamiltonian to 
simulate the effect of correlations. This term has the 
form 



SV(r) = -. 



2(r4 + r*) 



(4) 



where r$ is a cut-off parameter (we use ro = 1 as) and d 
is dipole polarizability of the core. We treat c as a fitting 
parameter and chose it to fit the known polarizability of 
dysprosium's ground state (166ag [24]]), which results in 
d = 3.7a%. 

Then we perform similar calculations for the excited 
states of the 4/ 9 6s 2 5d configuration, resulting in a cal- 
culated value of the static polarizability of 114 a^. Note 
that this approach does not distinguish between different 
states of the same configuration. Therefore, static polar- 
izabilities of all these states are assumed to be equal. This 
is only true for the static polarizabilities. Dynamic po- 
larizabilities are different for different states due to con- 
tributions of the near resonant states in Eq. ©. 

In the present paper we consider dynamic polarizabil- 
ities of four states of dysprosium: the even ground state 
(GS) and three odd long-lived excited states. The first 
excited state is 7 Hg at A = 1322 nm (£'=7565.60 cm" 1 ), 



TABLE I: Electric dipole transition amplitudes (reduced ma- 
trix elements in atomic units) used for calculation of the dy- 
namic polarizability of the Dy ground state Is- 



State n 




(cm- 1 ) 


(a.u.) 


4/ 9 5d6s 2 


'H 8 ° 


7565 


0.061 


4/ 9 5d6s 2 


7 H? 


8519 


0.124 


4/ 9 5d6s 2 


7to 

-M) 


9990 


0.059 


4/ 9 5d6s 2 


7,o 
J-8 


12007 


0.573 


4/ 9 5d6s 2 


7 G? 


12655 


0.108 


4/ 9 5d6s 2 


5 Kg 


13495 


0.424 


4/ 9 5d6s 2 


7to 
J-7 


14367 


0.475 


4/ 9 5d6s 2 


■ r ) T O 

J-8 


14625 


1.828 


4/ 9 5d6s 2 


5 H? 


15194 


1.452 


4/ 10 6s6p 


(8,0)1 


15567 


0.464 


4/ 10 6s6p 


(8,1)1 


15972 


1.365 


4/ 9 5d6s 2 


7 K1 


16288 


0.182 


4/ 10 6s6p 


(8,1)? 


16693 


1.842 


4/ 10 6s6p 


(8,1)1 


16733 


0.633 


4/ 9 5d6s 2 


7 K§ 


16717 


0.415 


4/ 9 5d6s 2 


7 K? 


17687 


0.763 


4/ 10 6s6p 


(8,2)| 


17727 


0.897 


4/ 10 6s6p 


(8,2)1 


18021 


0.684 


4/ 10 6s6p 


(8,2)? 


18433 


0.636 


4/ 9 5d 2 6s 


9 G? 


18528 


0.067 


4/ 9 5d 2 6s 


7 H 9 ° 


19557 


0.036 


4/ 9 5d6s 2 


5 K1 


19688 


0.627 


4/ 9 5d 2 6s 


7/-io 


21540 


0.523 


4/ 10 6s6p 


(7,2)| 


21838 


0.513 


4/ 9 5d 2 6s 


?° 


23271 


0.003 


4/ 10 6s6p 


(8,1)1 


23737 


12.277 


a A„ a = <n||D||o>. 



and we denote it as 01 for reference. This state is 
in the telecommunications band, and could be used for 
hybrid atom-photon telecom quantum information net- 
works. The second excited state is the 7 Ig state at 
1001 nm (9990.95 cm -1 ), we denote as 02. InAs quan- 
tum dots (QDs) emit in this wavelength range, allow- 
ing the possibility for hybrid quantum circuits of QD 
single photon emitters coupled to neutral atom-based 
long-lived quantum memory. 03 is the 5 Kg state at 741 
nm (13495.92 cm" 1 ), which is a closed cycling transition 
with a linewidth |28| optimal for creating a narrow-line 
MOT. States 02 and 03 could also be useful for resolved- 
sideband cooling, as discussed below. 

We calculate dynamic polarizabilities using the Eq. (fTJ) 
in which we substitute transition amplitudes found from 
the CI calculations [5j and experimental energies. We use 
theoretical values in the few cases where experimental en- 
ergies are not available. Tables HI and Irfl show calculated 
electric dipole transition amplitudes (reduced matrix el- 
ements) used in the calculations. The data from Table HT1 
can be used to calculate the lifetimes of the three excited 
states. The results are 5.2 ms for the 01, 2.7 ms for 02, 
and 21 /is for 03, though 03 has recently been measured 
to be 89.3 fis 
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TABLE II: Electric dipole transition amplitudes (reduced ma- 
trix elements in atomic units) used for calculation of the dy- 
namic polarizabilities of the selected three long-lived Dy ex- 
cited states. 



300 



State 


n 


-En 


\A na \ a (a.u 


) 






(cm- 1 ) 


01 b 


02 c 


03 d 


4/ lu 6s* 


is 





0.061 


0.059 


0.424 


4/ 10 6s 2 


5 l7 


4134 


0.007 






4/ 10 5d6s 


3 [8] 9 


17515 


0.033 


0.154 


0.139 


4/ 10 5d6s 


3 [7]s 


17613 


0.195 


0.321 


0.046 


4/ 10 5d6s 


3 [6] 7 


18095 


0.170 






4/ 10 5d6s 


3 [9] io 


18463 




0.060 


0.210 


4/ 10 5d6s 


3 [8]s 


18903 


0.300 


0.401 


0.268 


4/ 10 5d6s 


3 [7]r 


18938 


0.259 






4/ 10 6s 2 


3 K2 8 


19019 


0.083 


0.093 


0.064 


4/ 10 5rf6s 


3 [9]9 


19241 


0.079 


0.119 


0.015 


4/ 10 5d6s 


3 [10] io 
3 [9]s 


19798 




0.049 


0.192 


4/ 10 5d6s 


20194 


0.358 


0.336 


0.288 


4/ 10 5d6s 


3 [10] 9 

/IS Is 
{-£> 2> 7 


20209 


0.041 


0.062 


0.010 


4/ 9 6s 2 6p 


20614 


2.164 






4/ 9 6s 2 6p 


( 15 1 \ 


20790 


3.737 


4.254 


2.911 


4/ 10 5d6s 


21074 


0.550 






4/ 10 5d6s 


3 [7]s 


21603 


0.362 


0.450 


0.130 


4/ 10 5d6s 


3 [6]r 


21778 


0.463 






4/ 10 5d6s 


1 [9] 9 


22046 


0.079 


0.081 


0.116 


4/ 10 5rf6s 




22487 




0.004 


0.112 


4/ 9 5d6s6p 


?9 


23219 


1.313 


3.879 


0.400 


4/ 10 5d6s 


?7 


23361 


0.172 







a A na = (n||D||o); 
'State a = 4/ 9 5d6s 2 , 
c State a = 4f 9 5d6s 2 , 
d State a = 4/ 9 5d6s 2 , 



7tto 

H g , 

7to 
L 9< 

5 K° 



E - 

E -- 

E : 



: 7565.60 cm" 
9990.95 cm~ 
: 13495.92 cm 



A = 1322 nm; 
A = 1001 nm; 
A = 741 nm. 



Figures [TJ [5] and [3] show dynamic polarizabilities of 
three pairs of states: GS and 01 (Fig. [T]), GS and 02 
(Fig. [2]) and GS and 03 (Fig. [3]). Lines crossing means 
that energy shifts of two states in a laser field are identi- 
cal, and atoms have the same oscillation frequency in op- 
tical dipole traps at this wavelength regardless of whether 
they are in their ground or excited state. These so-called 
magic wavelengths occur most often very close to narrow 
resonances. 



B. Simple estimations 

In this subsection we present a way of estimating magic 
wavelengths for complex atoms in the vicinity of narrow 
resonances. Although, all magic wavelengths presented 
in this work are found by the many-body calculations, 
the formulas of this subsection can be used to find more 
magic wavelengths for dysprosium or to find magic wave- 
lengths for other complex atoms. We demonstrate that 
many-body calculations and simple estimations give close 
results. 

In the case of a narrow resonance, energy denominator 
in Eq. ([!} is close to zero. This makes it possible to write 
an approximate formula for the magic frequency corre- 



200 



100 



12.6 



12.8 13 13.2 

a (10 3 cm -1 ) 



13.4 



FIG. 1: Dynamic polarizability a of the ground state of Dy 
(solid line) and 01 (dotted line) between laser frequencies 
(wavelengths) 12500 cm" 1 (800 nm) and 13500 cm -1 (741 
nm). Lines cross at magic frequencies, and large dots corre- 
spond to the most useful. 



300 



200 




100 - 



10 11 12 

cj (10 3 cm" 1 ) 



13 



FIG. 2: Dynamic polarizability a of the ground state of Dy 
(solid line) and 02 (dotted line) between laser frequencies 
(wavelengths) 8500 cm" 1 (1.18 ^m) and 13000 cm" 1 (769 
nm). Lines cross at magic frequencies, and large dots corre- 
spond to the most useful. 



sponding to this wavelength. Starting from the condition 
acs(w*) = a a (uj*), (5) 
and presenting dynamic polarizability of the excited state 
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300 



200 - 



100 




FIG. 3: Dynamic polarizability a of the ground state of Dy 
(solid line) and 03 (dotted line). Lines cross at magic fre- 
quencies, and the large dot corresponds to the most useful. 



in the vicinity of the resonance n in the form 



a a (ui) = a a (0) 



1 



3(2J a + 1) 



1 



E a — E n 



1 



a||D||n} 2 , 



E a — E n — uj 
we arrive at the following expression: 

E n — E 



E. n 



Sn = 



1 



\E a — E n \ 
(a||D||n) 2 



(6) 



(7) 



3(2J Q + l)(a GS (w„)-a Q (0)) : 



where w„ = \E a — E n | is the resonance frequency, n is the 
resonance number, (<z||D||n) is the electric dipole transi- 
tion amplitude from the excited state a to the resonance 
state n, and acs(u n ) is the scalar polarizability of the 
ground state at w„. Note that uj n is the resonance fre- 
quency for the upper state, and the polarizability of the 
ground state usually changes very little in the vicinity of 
the uj n 's. 

In case of two closely spaced resonances in the upper 
state polarizability, magic frequencies in the vicinity of 
resonance energies Ei and E 2 can be found using approx- 
imate formula 



where 



S12 = A 



12" 



a 



Ei + E2 



05(^12) 



<5l2 



„C°) 



(0) 



A\ 



A\ 



A\ 



A\ 



(8) 



(9) 



Here, A 12 = {E x - E 2 )/2, Cj = 3(2J Q + 1), uj 12 = (E 1 + 
E 2 )/2, A x = (a||D||l) and A 2 = (o||D||2). 

One does not need to know dynamic polarizability of 
an excited state to find magic frequencies using Eqs. Q 
or (JSJ). However, one still needs to know the dynamic 
polarizability of the ground state and the relevant tran- 
sition amplitudes. An approximate solution can be found 
using the following rules: 

• The dynamic polarizability of the ground state is 
fitted very well within energy interval from to 
0.05 a.u. by 



a GS (uj) 



2.955 



2.955 



u+0. 10815 u-0. 10815 

110 + 3000w 2 . 



(10) 



We keep here numerical parameters for the domi- 
nant contribution to aosi which is due to the tran- 
sition to the 421-nm state with energy E= 23737 
cm -1 = 0.10815 a.u. This transition dominates 
due to the largest value of the transition amplitude: 
(a||D||n) =12.277 a.u. (see Table H|). The last two 
terms in Eq. (fTU)) fit the contribution of all other 
transitions in the energy interval < uj < 0.05 a.u. 
Resonances within this interval are ignored since 
they are too narrow for any practical importance. 

• Transition amplitudes can be estimated using ap- 
proximate selection rules. If the difference of the 
total angular momentum L between two states is 
larger than 1, or if total spin S of the states is not 
equal, then the amplitude is not zero due to rela- 
tivistic corrections but is likely to be small. One 
can use for a rough estimate A=0.1 a.u. A simi- 
lar estimation can be used if the transition is sup- 
pressed by configuration mixing, i.e. the transition 
between leading configurations cannot be reduced 
to a single-electron allowed electric dipole transi- 
tion. If no selection rules are broken, the amplitude 
is likely to be large and one can use A=3 a.u. as a 
rough estimate. 

• Polarizability of the excited state at zero frequency 
can be estimated using Eq. (fTJ) with experimental 
energies and with the amplitude estimated using 
the procedure in the previous paragraph. 

This procedure can help in estimating magic wavelengths 
not only for the transitions considered in present paper 
but also for some other transitions. The main condition 
for it to work is that magic wavelength should be close to 
a resonance so that resonance term dominates in Eq. ([1} . 



III. RESULTS 

Table HIT] shows the magic wavelengths and correspond- 
ing polarizabilities for the three transitions in Dy. We 
substitute calculated transition amplitudes from Table [TT1 
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TABLE III: Magic wavelengths (A*), frequencies (u)*, cm 1 ), and polarizabilities (a) for the three transitions in Dy. 



Transition 


Resonance 




Magic frequencies 




A* 


a 




E n — E a 


<>n 


S n c Formula 


Calculations 


(nm) 


(a.u.) 


GS - 01 


11337 


4 


11333 


11326 


883 


192 


£ a =7566 cm" 1 


11372 


3 


11368 


11366 


880 


192 


A = 1322 nm 


13224 


537 


12686 


12638 


791 


202 




13136 




-35 13101 


13100 


763 


208 


GS - 02 


7523 


1 


7522 


7521 


1330 


172 


£ a =9991 cm" 1 


7622 


4 


7617 


7613 


1314 


174 


A = 1001 nm 


8912 


7 


8905 


8891 


1125 


178 










9749 


1026 


181 




1 Q007 

16ZZ ( 


/ICQ 


1 07/1 /I 


1 rin-i 7 

1/81 / 


/ oU 




GS - 03 


5407 


4 


5403 


5401 


1850 


163 


£ a =13496 cm" 1 


6302 


2 


6300 


6297 


1588 


171 


A = 741 nm 


6697 


4 


6693 


6671 


1494 


172 




7293 


479 


6814 


6812 


1468 


172 




9722 


8 


9714 


9716 


1029 


181 



"Using formula 10, ui* n = E n — E a — 8 n ; 
b E n = (20614 + 20790)/2 - 7566; 

c Using formulas J§E), oj* 12 = (E 1 + E 2 )/2 - E a + 5 12 . 



when using the analytical formulas (0 and ((8|). The col- 
umn marked as calculations presents magic frequencies 
which come from numerical calculations. Magic wave- 
lengths in the next column correspond to calculated fre- 
quencies. Most of the magic frequencies are due to very 
narrow resonances and might be inconvenient for prac- 
tical use due to optical dipole trap frequency instabili- 
ties and enhanced spontaneous emission. However, there 
are magic wavelengths for each of the three transitions 
where the resonance is not very narrow or even absent. 
They are A = 791 nm (uj* = 12638 cm" 1 ) and A = 763 
nm (uj* = 13100 cm" 1 ) for the 1322-nm GS-01 transi- 
tion; A = 1125 nm (uj* = 8891 cm" 1 ), A = 1026 nm 
(uj* = 9749 cm" 1 ) and A = 780 nm (uj* = 12817 cm" 1 ) 
for the 1001-nm GS-02 transition; and A = 1029 nm 
(uj* = 9716 cm" 1 ) and A = 1468 nm (uj* = 6812 cm" 1 ) 
for the 741-nm GS-03 transition. 

The magic frequency uj* = 13100 cm" 1 for the GS- 
01 transition is between two close resonances at uj = 
13149 cm" 1 and uj* = 13224 cm" 1 which correspond 
to transitions from the 7 H§ state at E = 7566 cm 
to the close states of the 4f 9 6s 2 6p 1 / 2 configuration at 
E = 20614 cm" 1 and E = 20789 cm" 1 . Using formu- 
las (|8]|9|) gives a very accurate estimate of the magic fre- 
quency (see Table IHl]) . It is interesting to note that there 
is a frequency interval for the GS-01 transition where po- 
larizabilities of two states come very close to each other 
but don't cross: Sa/a < 2% for 12121 < uj < 12183 
cm" 1 . 

The magic wavelengths A = 780 and 1026 nm for the 
GS-02 transition and A = 1468 nm for GS-03 do not 
correspond to any very-close resonance, and the values 
of the polarizabilities coincide by chance rather than due 
to a resonance. Therefore, these magic wavelengths are 



the least sensitive to laser frequency fluctuations and are 
most promising for resolved-sideband cooling, precision 
measurement, and QIP applications. The GS-02 tran- 
sition magic wavelengths can be reached with high opti- 
cal power using a Tksapphire or tapered amplified diode 
laser for 780 nm and diode laser or fiber laser for 1026 
nm. The GS-03 transition magic wavelength at 1468 
nm could be reached with a low-power diode laser, which 
could perhaps be doped-fiber amplified, and the 1029-nm 
wavelength could be reached with a fiber laser. 

Future work will explore the role of laser polarization 
on the magic wavelength position. Such a calculation is 
beyond the scope of this present work, but we estimate 
that the shifts will be small since the magic wavelengths 
for unpolarized light occur in proximity to resonances. 

Optical dipole traps at these wavelengths, far from the 
broad Dy transitions, would be suitable for lattice con- 
finement in the Lamb-Dicke regime without undue heat- 
ing. For example, ID lattice confinement at the 780-nm 
magic wavelength with 0.5 W provides ample trap depth 
with sub-1 Hz scattering rates. With larger laser inten- 
sities, suitable trap depths and low scattering rates can 
be achieved at the other magic wavelengths. Vibrational 
spacing can be many tens of kHz, which is large enough 
for resolved-sideband cooling on the 2 kHz-wide 741-nm 
transition [28|. For rapid cooling, the 50 Hz- wide 1001- 
nm transition — much narrower than any trap frequen- 
cies in a typical 3D optical lattice — would need to be 
broadened via a quenching transition [22j, and the opti- 
cal dipole trap magic wavelength would need to be ad- 
justed to compensate the Stark shift from the quenching 
laser. Resolved-sideband cooling on these narrow tran- 
sitions in a 3D optical lattice may provide an alterna- 
tive route to quantum degeneracy [29( versus evapora- 
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tive cooling, which may fail due to (as yet unmeasured) 
unfavorable scattering properties in this highly dipolar 
gas. 
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